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摘 要 针对 工程 中 的 反问 题 ,以 POD (最 佳 正 交 分 解 ，Proper Orthogonal Decomposition) 技术 为 基础 设计 了 一 种 
反问 题 的 低 阶 计算 模型 ， 通 过 采集 到 的 有 限 输出 结果 ， 仅 采用 简单 的 非 线 性 规划 与 插值 法 即 可 求解 该 模型 ， 大 幅度 
地 降低 了 工程 中 的 反问 题 的 求解 难度 与 计算 时 间 。 由 于 模型 求解 过 程 不 涉及 原 物理 问题 的 控制 方程 ， 因 此 ， 该 模型 
在 实际 中 具有 一 定 通用 性 。 以 圆 管 内 强制 对 流 换 热 问题 为 应 用 实例 ， 结 果 表 明 ， 基 于 一 些 热电 偶 测 量 的 温度 数据 ， 
低 阶 反问 题 计算 模型 能 够 较 精确 地 反 算 出 未 知 的 流体 进口 温度 与 圆 管 外 壁 热流 密度 , 与 真实 值 相 比 仅 相差 1.0% 左 右 。 
此 外 ， 与 经 典 的 共 轿 梯度 法 相 比 ， 低 阶 反问 题 计算 模型 的 计算 速度 可 提高 1200 倍 以 上 。 
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Design of reduced-order computational model of reverse problem using POD 
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Abstract This paper proposes a reduced-order computational model for reverse problem mainly based on 
POD (Proper Orthogonal Decomposition), through finite output data collected in advance, the model is 
solved only using non-linear programming and interpolation methods, which considerably reduces the 
difficulty and computational time of reverse problems in engineering. The governing equations of the 
original problem are not considered in the model, which suggests the universality of the reduced-order 
reverse problem computational model. The reduced-order reverse problem computational model is 
approved on a forced convection and heat transfer problem in a circular pipe, the results show that the 
reduced-order reserve problem computational model enables us to accurately calculate the unknown inlet 
temperature and wall heat flux based on some temperature data measured by thermocouples, where the 
error between the computational results and the actual results is only about 1.0%. Besides, it is quite fast, 
where the achieved increase in calculation speed is more than 1200 times compared with the classical 
conjugate gradient methodology. 
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工程 中 的 反问 题 ， 或 称 逆 问题 ， 其 基本 思想 是 
根据 输出 结果 的 部 分 信息 反 算 出 模型 的 未 知 参数 。 
反问 题 的 出 现 也 是 实际 工程 问题 的 需要 ， 特 别 是 许 
多 由 于 条 件 限制 而 不 易 对 想 要 知道 的 参数 进行 直接 
测量 的 情况 ， 如 预测 未 知 材料 的 物性 、 油 藏 中 岩石 
的 物性 、 结 构 工 程 中 的 损坏 探测 、 传 热 问 题 中 边界 
温度 或 热流 密度 等 参数 的 监测 等 。 反 问题 的 求解 方 
法 中 ， 以 共 斩 梯度 法 口 应 用 得 最 为 广泛 ， 共 罗 梯 度 法 
要 求 计算 上 且 标 消 数 对 所 求 参数 的 导数 《 称 为 敏感 系 
数 )， 再 借助 于 敏感 系数 而 逐渐 台 近 所 要 求解 的 模型 
参数 。 由 于 采用 共 斩 梯 度 法 求解 敏感 系数 需要 借助 
于 及 反复 求解 原 问题 的 控制 方程 ， 这 需要 消耗 大 量 
的 计算 时 间 与 存储 空间 ， 从 而 使 得 它 无 法 应 用 于 那 
些 需 要 实时 预测 未 知 参数 的 领域 。 此 外 ， 工 程 中 许 
多 问题 的 物理 过 程 由 一 组 极其 复杂 的 非 线 性 偏 微 分 
方程 所 控制 ， 很 难 甚至 无 法 求解 敏感 系数 。 

POD 〈 最 佳 正 交 分 解 ， Proper Orthogonal 
Decomposition) “技术 是 从 已 知 的 对 应 于 若干 个 参 
数 〈 控 制 参数 或 物性 参数 ) 的 物理 场 样本 (由 数值 
模拟 或 实验 获得 ) 中 提取 出 一 组 只 与 空间 变量 有 关 
的 正 交 基 函 数 ， 这 些 基 函数 代表 了 原 物理 系统 的 动 
力学 特性 ， 再 与 Galerkin “投影 方法 结合 即 可 实现 对 
原 物 理 系 统 控制 方程 的 降 阶 ， 即 建立 低 阶 模型 ， 这 
里 的 “ 阶 ”代表 每 个 时 间 步 长 内 所 要 求解 变量 的 个 
数 ， 利 用 这 种 低 阶 模型 可 快速 地 计算 出 参数 可 变化 
范围 内 的 任 一 参数 所 对 应 的 物理 场 。 这 也 是 目前 以 
POD 作为 降 阶 工具 的 主流 降 阶 方法 ， 也 是 其 最 典型 
的 应 用 。 通 过 该 方法 ，POD 已 被 成 功 应 用 于 许多 科 
学 与 工程 领域 ， 如 信号 分 析 趾 、 流 体 流动 与 传 热 、 
结构 动力 学 说 、 最 优 控制 四 以 及 反问 题 呈 六 等 。 对 
于 反问 题 ， 文 献 [12] 与 [13] 分 别 将 低 阶 模型 引入 圆 管 
内 层 流 对 流 换 热 及 饱和 地 下 水 流动 的 反问 题 模型 ， 
求解 方法 的 主体 依然 是 共 轿 梯度 法 ，POD 在 每 一 次 
求解 敏感 系数 过 程 中 对 所 需求 解 的 控制 方程 实现 降 
阶 ， 从 而 降低 计算 时 间 ， 尽 管 最 终 计 算 时 间 降 低 了 
数 十 倍 ， 但 这 仍然 属于 主流 降 阶 范畴 ，POD 的 降低 
计算 时 间作 用 仅 在 求解 控制 方程 这 一 步 得 到 发 挥 ， 
且 实施 难度 随 着 控制 方程 复杂 度 的 增加 而 增加 ， 如 
文献 [12] 中 的 低 阶 模型 仅 适 用 于 层 流 ， 对 于 汕 流 便 无 
法 实施 。 
鉴于 上 述 问 题 ， 本 文采 用 POD 的 另 一 种 使 用 方 
法 , 即 POD 插值 法 设计 了 一 种 新 的 反问 题 计算 模型 ， 


Tl 


即 低 阶 反 问题 计算 模型 。 以 圆 管内 强制 对 流 换 热 为 
应 用 实例 对 该 方法 的 精确 性 与 计算 效率 进行 了 检 


$, 
1 POD 技术 简介 


以 w, (n=1 N) 代表 参数 向 量 , PN 为 向 
量 的 个 数 ，q(X,w, ) 代表 在 参数 yy 下 达到 稳 态 时 刻 
的 某 个 状态 参数 样本 , 其 中 x 代表 空间 变量 。 从 分 别 
对 应 于 参数 w, (n =1… N) N 个 状态 参数 样本 中 
提取 出 一 组 只 与 空间 变量 有 关 的 正 交 基 函数 
办 C) (m 2 1-- M) 后， 状态 参数 g(x,W,) 可 通过 式 
(D 得 到 重 构 ， 即 : 


M 
q(x, w,)=} a, w, p(x n=LN D 
m-l 
式 中 ，M 为 截断 自由 度 ， 即 所 采用 的 基 函 数 的 数量 ; 
Cs Ow, ) 为 对 应 于 各 基 函 数 的 系数 ，C, (Ww ) 的 大 小 
RETE m^ POD ERBA, O) 对 重 构 g(x,w,) 所 
贡献 的 程度 ， 对 于 稳 态 时 刻 ， 只 与 参数 w, 相关 。 
1.1 “快照 ”方法 
POD 基 函 数 的 求解 可 由 Sirovich" ^32 t fre 
照 "方法 完成 。Sirovich 将 POD 基 函 数 表 示 为 各 样本 
(TES, B. 


N 
6, (x) 2 3 P q(xw,)m-l-M D 
n-l 
RP, RAD 代表 第 n 个 参数 w, 对 应 的 样本 对 构 
建 第 m 个 POD 基 函 数 的 贡献 程度 , 可 由 求解 如 下 特 
征 值 问 题 获得 : 


Cb" =A"b" (3) 


式 中 ，C 为 一 个 N 维 的 对 称 矩 阵 ; A" JC 的 第 n 个 
特征 值 , b” 为 对 应 的 特征 向 量 , p" 的 元 素 便 为 式 (4) 


中 的 系数 bp" o JEREC 的 元 素 为 : 
c =}f Qr, w,)q(x,w dQ 
hi N af YU qxx, j (4) 


1.2 能 量 最 优 性 
X G) 中 的 特征 值 4” 大 小 代表 了 相应 的 POD 
函数 在 原 物 理 场 中 所 占据 的 能 量 多 少 。 定 义 参 数 : 


b 


M N 
$20» 4 (5) 
n=l n=l 
WERI, 23 5, 表示 前 M 个 POD 基 函 数 对 于 总 
体能 量 的 贡献 程度 。 
根据 POD 基 函 数 的 “能 量 最 优 性 ”3 原则 ， 按 
特征 值 的 数值 由 大 到 小 排列 ， 仅 前 M (MN) 个 特 
征 值 便 能 使 参数 6 接近 1。 正 是 由 于 POD 基 函 数 的 
这 种 能 量 最 优 性 , 使 重 构 公 式 (1) 仅 使 用 极 少 的 POD 
基 函 数 便 可 较 高 精度 地 重 构 出 原状 态 参数 。 
1.3 POD 插值 
对 于 设计 参数 区 , 下 的 系数 w Ov.) 可 由 如 下 的 
投影 公式 得 到 ， 即 : 


a, Qw,) = (6, (7),4(x,w,)) 
m-l--M,nzl-N 


对 于 非 设计 参数 ww 下 的 系数 Q,,(w) 可 由 设计 参 
Aw, 下 的 系数 Q,, (w,) 进行 插值 Poe, BI 


QW)=f (0 Ww)) m= M,n=1.N (7) 
式 中 ， 函 数 了 由 所 采用 的 插值 方式 决定 。 
2 低 阶 反问 题 计算 模型 


2.1 数学 模型 

基于 POD 插值 方法 ， 本 文 所 提出 的 低 阶 反问 题 
计算 模型 的 基本 思想 为 : 根据 稳 态 时 刻 测量 到 的 物 
理 场 信息 q ， 求 解 出 其 对 应 的 POD 系数 
a (m=1…M) ， 再 从 一 系列 设计 参数 
w, (nl N) 中 采用 插值 法 ( 式 CDD 寻找 出 w 
所 对 应 的 参数 w* 。 其 数学 模型 可 写 为 : 


(6) 


min 


(8) 


lo = > Anpa (x) 


m=1 


2 


st. cn <a, xa"* (m=1 M) (9) 


min Se -a. 00] (10) 


m= 


= 


s.t. a, (w) = f (a, (w,)) (nzLl-N) (10D 


RP, Iure Rt 2- 范 数 。 

R (8) 与 式 (9) 用 于 计算 出 对 应 于 稳 态 时 刻 
已 知 物理 场 9 H POD 3c o7, (m=1…M)， 其 中 
Han, tg" 分 别 为 POD 系数 的 取 值 下 限 与 上 


限 ， 由 一 系列 设计 参数 ww, (n =1…N) 对 应 的 POD 
系数 取 值 范围 决定 。 不 难 发 现 ， 求 解 
a, (m - 1 M) 的 过 程 只 是 个 简单 的 非 线 性 规划 
问题 。 式 (10) 与 式 (11) 为 采用 插值 法 从 一 系列 设 
计 参 数 w, (n=1…N) rms dd a, (m=1…M) 
所 对 应 的 最 优 参 数 w 的 计算 模型 ， 计 算 中 ， 只 要 
Qt, Ow) POSER SEU 使 得 式 (10) 满足 一 定 的 求 
解 精度 ， 则 此 时 的 参数 便 为 所 要 寻找 的 最 优 参数 


W o 
22 计算 步骤 

本 文 设计 的 基于 POD 技术 的 低 阶 反问 题 计算 模 
型 的 计算 步骤 可 总 结 如 下 : 

COD 在 参数 取 值 范围 内 均匀 选取 N 个 设计 参数 
f& w, (n=1…N)，, 计算 每 个 参数 值 对 应 于 稳 态 时 刻 
某 些 测量 点 的 状态 参数 q(x,w,) n1 N), 3k 
得 NN 个 样本 ; 

(2) 对 样本 实施 POD 技术 求 取 一 组 正 交 基 函 
9, (x) (m=1…M); 

(3) 由 式 (6) 计算 出 所 有 设计 参数 下 的 POD 
系数 (w,) (m=1…M,n=1…N)， 并 确定 最 小 
值 a™ 5E" s 

(4) 利用 式 (8) 与 式 (9) 的 非 线 性 规划 模型 
计算 出 对 应 于 4 和 的 最 优 POD 系数 
a^, (mi =1 M); 

(5) 从 一 系列 设计 参数 w, (n=1- N) 中 利用 
X 1) 的 插值 公式 寻找 出 使 得 式 (10) 满足 一 定 精 
度 的 最 优 参 数 w* ， 其 具体 的 求解 步骤 为 : 

Q 为 式 (10) 设置 求解 精度 &; 

© 在 计算 区 域内 将 参数 的 取 值 数量 加 倍 为 
w, (nzl--2N); 

© AAR aD 逐个 计算 出 新 增 的 参数 值 所 对 
应 的 POD 系数 ; 

QD 将 各 个 参数 对 应 的 POD 系数 逐个 代入 式 
(100, dX (100 满足 精度 & ， 则 计算 结果 收敛 ， 
此 时 的 POD 系数 (w) (m =1…M) 所 对 应 的 参数 
便 是 所 要 求解 的 最 优 参 数 w* ， 若 不 满足 ， 则 由 步 双 
@ 开 始 重新 计算 ， 直 至 式 OO 满足 精度 & 。 


需要 指出 的 是 , 在 pa -a,(w)] 极 小 值 
m-l 


点 的 某 个 邻 域内 可 能 存在 多 个 参数 所 对 应 的 POD 系 
数 满足 精度 & , 这 时 只 需 再 从 中 选取 使 得 目标 函数 取 
得 最 小 值 的 参数 即 可 。 

由 上 述 求解 步骤 可 看 出 ， 与 文献 [12] 与 [13] 中 的 


POD 主流 降 阶 方法 相 比 ， 低 阶 反 问题 计算 模型 不 涉 
及 控制 方程 , 而 只 以 少量 的 POD 系数 作为 中 间 变 量 ， 
仅 采 用 简单 的 非 线性 规划 与 插值 方法 即 可 求解 出 所 
要 寻找 的 参数 ， 实 现 了 直接 对 反问 题 进行 降 阶 ， 这 
可 大 幅度 地 降低 反问 题 的 求解 难度 与 计算 时 间 。 此 
外 ， 由 于 模型 求解 过 程 不 涉及 控制 方程 ， 即 和 所 研 
究 问 题 复杂 与 否 无 关 ， 因 此 ， 本 文 所 提出 的 新 方法 
在 实际 工程 应 用 中 具有 一 定 的 通用 性 。 另 外 ， 值 得 
指出 的 是 ， 本 文 提出 的 低 阶 反问 题 计算 模型 只 适用 
于 求解 待 求解 参数 为 定 值 的 问题 ， 即 不 随时 间或 空 
间 而 变化 。 


3 ”应 用 实例 


3.1 问题 描述 

本 文 以 文献 [12] 中 的 圆 管内 稳 态 强制 对 流 换 热 
问题 为 应 用 实例 具体 说 明 低 阶 反问 题 计 算 模 型 的 实 
施 过 程 ， 其 模型 如 图 1 所 示 。 流 体 入 口 温 度 T, KA 
管 外 壁 热流 密度 9 为 待 求解 参数 ， 且 为 定 值 ; 圆 管 半 
hr, 20.01 m, 长 度 为 1=1m , 不计 管 壁 导 热 热 阻 ; 
流体 导热 系数 为 n=0.599W/(m:K) ， 密 度 为 


ME 


IIi] 


p=1000kg/ m, ERAH c, =4.183 kJ/(kg: K), 
流体 入 口 流速 沿 径 向 的 分 布 为 : 
u(r) -2u,0 - 7? /r?) (12) 


AP, US 为 平均 流速 ，0.05 m/s. 

同 管 内 的 黑 点 代表 热电 偶 ， 假 设 沿 轴 癌 均 匀 地 
布置 了 25 个 热电 偶 用 于 测量 稳 态 时 的 流体 温度 ， 
为 测量 点 位 置 与 壁面 的 径 向 距离 ， 不 计 测 量 误差 。 
现 基于 热电 侦 测 量 的 温度 场 采用 本 文 提出 的 低 阶 反 
问题 计算 模型 反 算出 流体 入 口 温 度 7, 及 圆 管 外 壁 热 
流 密度 49 。 不 计 重力 影响 ， 模 型 网 格 划分 为 2000x 
100, 


这 里 需要 说 明 的 是 ， 尽 管 本 文 提 出 的 低 阶 反问 
题 计 算 模 型 求解 过 程 不 涉及 控制 方程 ， 即 对 于 本 文 
的 应 用 实例 无 论 是 层 流 还 是 满 流 均 适用 ， 但 之 所 以 
本 文采 用 层 流 ， 主 要 是 为 了 对 比 本 文 新 方法 与 经 
的 共 轿 梯度 法 计算 效率 的 需要 ,文献 [15] 已 给 出 采 月 
共 轿 梯度 法 求解 本 文 应 用 实例 的 具体 实施 过 程 。 而 
对 于 满 流 ， 由 于 控制 方程 的 极其 复杂 性 ， 共 轿 梯度 
法 几乎 无 法 实施 ， 目 前 也 并 无 相关 文献 报道 。 
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1 物理 模型 


Fig. 1 Physical model 


3.2 POD ŽAK 

ER 2 个 待 求解 参数 的 取 值 范围 及 取 值 情况 如 
表 1 所 示 ， 将 各 个 参数 值 任意 搭配 共 得 到 48 个 参数 
组 合 。 采 集 每 个 参数 组 合 下 稳 态 时 在 轴线 


rfr, 20.2. rfr, 20.15 & r. /r=0.1 上 热电 偶 所 
在 网 格 处 的 温度 值 ， 共 得 到 48 个 样本 。 对 样本 实施 
POD 技术 ， 表 2 给 出 了 前 5 个 特征 值 及 其 相应 的 能 
量 分 布 。 


表 1 待 求 参数 的 取 值 范围 及 取 值 情况 


Table 1 The ranges and values of parameters to be Solved 


空 制 参数 取 值 范围 取 值 
入 口 温度 288-313 K 288, 293, 298, 303, 308, 313 


热流 密度 3000-10000 W/m? 


3000, 4000, 5000, 6000, 7000, 8000, 9000, 10000 


3k 2 前 5 个 特征 值 及 对 应 的 能 量 分 布 


Table 2 The first five eigenvalues and corresponding energy distribution 


n 1 2 3 4 5 
A 23978.9272 80.6821 3.2631 0. 2162 0.0136 
Č u (9) 99.6425 99.9778 99.9914 99.9923 99.9924 
由 表 2 可 见 仅 前 5 个 基 函 数 便 几 乎 捕捉 到 了 样 Fig. 2 Reconstruction error for partial samples 
本 的 全 部 能 量 份额 。 3.3 POD ZEER IUE GR E 1r 
图 2 为 部 分 样本 的 重 构 误 差 。 由 图 2 可 以 看 出 ， 图 3 与 图 4 分 别 给 出 了 由 投影 公式 ( 式 (6)) 
重 构 误差 随 着 截断 自由 度 M. 的 增加 而 迅速 减 小 ， 计算 出 的 各 样本 的 第 1、2 组 POD 系数 ， 其 中 ， 图 
并 趋 近 于 一 个 定 值 。 当 M=5 时 ,五 的 值 已 降 至 0.05% (a) 中 的 每 个 变化 周期 内 的 样本 是 按 入 口 温 度 不 变 ， 
左右 。 其 中 五 的 定义 为 ; 外 壁 热流 密度 不 断 增加 的 方式 安排 ， 图 Cb) 则 是 外 
la -a,,], 壁 热流 密度 不 变 ， 而 入 口 温 度 不 断 增 加 。 可 以 看 出 ， 
E =——  x10076 (13) ”样本 POD 系数 随 着 参数 的 变化 几乎 呈 线 性 变化 ， 其 
lal; È 3 组 POD 系数 同样 存在 类 似 的 变化 关系 ， 这 里 不 


NF, d 为 数值 求解 结果 ，9x 为 重 构 结果 。 给 出 。 因 此 ， 本 文采 用 线性 插值 方法 计算 非 
SMS 设计 参数 所 对 应 的 POD 系数 。 
外 壁 热流 密度 : 
—ma— 293 K, 4000 W/m? 
—o— 303 K, 6000 W/m? 
—4— 308 K, 8000 W/m? 
04 
$ 03 
y 
0.2 
0.1 
全 一 
0.0 T T T T T T T T T T 
1 2 3 4 5 6 7 8 9 10 
M 
图 2 部 分 样本 的 重 构 误差 
6004 16004 
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/ | d 
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4504 1450 - 
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| I 
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图 3 样本 的 第 1 组 POD 系数 
Fig. 3 The first group of POD coefficients for the collected samples 
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图 4 样本 的 第 2 组 POD 系数 
Fig. 4 The second group of POD coefficients for the collected samples 
为 了 检验 POD 基 函 数 用 于 非 设 计 参 数 下 插值 结 ” 流体 入 口 温度 为 305 K. 
果 的 精确 度 ， 图 5 给 出 了 部 分 非 设计 参数 下 在 轴线 由 图 5 可 见 ， 插 值 结果 与 数值 求解 结果 基本 
/n=0.2 上 的 线性 插值 结果 与 直接 数值 求解 结 致 , 对 于 图 5(a), 平均 误差 已 为 0.367%, 对 于 图 5(b)， 
的 对 比 ， 其 中 ， 选 取 了 前 5 个 POD 基 函 数 。 图 中 的 EX 0.413%, 说 明 所 求 得 的 POD 基 函 数 跨 越 了 原 问 


线条 为 数值 求解 结果 ， 几 何 图 标 为 插值 求解 结果 。 


J2s ix EG 2 en 
对 于 图 5(a)， 热 流 密度 为 7500 W/m ， 对 于 图 5(b)， 
3357 流体 入 口 温 度 
3304 i 290.5 K 
1 ae o 2905K 
3254 ge^ Le * o 205.5 K 

J 一 __ j 

3203 _ pi a< Y 2955K 
ga” pr NT. 300.5 K 
3154 p SEM EP 4 v—Y «4 3005K 


测量 点 


(a) 


题 参数 较 大 范围 内 的 解 所 处 的 空间 ， 可 在 参数 允许 
变化 范围 内 通过 插值 不 断 重 复 使 用 。 


图 5 非 设计 参数 下 的 插值 结果 与 数值 求解 结果 的 对 比 


Fig. 5 Comparison of interpolation and numerical results for off-design parameters 


3.4 计算 结果 

图 6 至 图 8 分 别 代表 热电 偶 在 轴线 六 /=0.2、 
/n=0.15 及 x/r,=0.1 上 测量 到 的 几 组 温度 妆 
Jn. de 3 给 出 了 低 阶 反问 题 计算 模 型 的 计算 结果 与 
真实 结果 的 对 比 情况 ， 其 中 ， 求 解 精度 为 =107 。 
表 3 可 以 看 出 所 求解 结果 与 真实 值 较 接 近 ， 平 均 
相对 误差 e 仅 为 1.069%， 其 中 ， 平 均 相 对 误差 e 的 
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] 外 壁 热流 密度 : 
330] X 4500 W/m? 
E o 4500 W/m? 
325 o 5500 W/m* 
] v 5500 W/m? 
3204 = 7500 W/m? 
T ] « 7500 W/m 
" 315] — 8500 W/m 
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NL 4g N 
i 公式 为 : 
lo. -o| 
= 一 (14) 
lel, 


AP, 0o, 为 真实 结果 ，w 为 计算 结果 。 


一 e 一 数据 1 
一 < 一 数据 3 
—e— 数据 4 
一 < 一 数据 5 


DO T a a a a a 
0 10 15 20 25 30 


测量 点 


图 6 r,/r, - 0.2 处 热电 偶 测 量 到 的 几 组 温度 数据 


Fig. 6 Sets of temperature data measured by thermocouples at the location of ra / r,-02 
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345 o— 数据 1 3504 数据 1 
340 一 vy 数据 2 3453 一 7 一 数据 2 
335 一 < 一 数据 3 3403 一 4 一 数据 3 
一 e 一 数据 4 3353 一 e 一 数据 4 
330 x 数据 E 一 x 一 数据 5 
325 due 
3254 
gon T 320} 
x 315 = 3154 
d i E 
zi 310 T 3104 
305 3054 
300 3004 
295 395] 
290 2904 
285 -rr mer T 一 一 一 T 285 4—————————————————————————————-4 
0 5 10 15 20 25 30 0 5 10 15 20 25 30 
测量 点 测量 点 
图 7 rn, /r, =0.15 处 热电 偶 测量 到 的 几 组 温度 数据 图 8 r,/r, = 0.1 处 热电 偶 测量 到 的 几 组 温度 数据 
Fig. 7 Sets oftemperature data measured by thermocouple Fig. 8 Sets oftemperature data measured by thermocouple 
at the location of 7, /r, — 0.15 at the location of 7/7, — 0.1 
表 3 低 阶 反问 题 计 算 模 型 的 计算 结果 与 真实 结果 的 对 比 
Table 3 Contrast between the results of inverse problem model and real results 
温度 场 数 计算 结果 真实 结果 a - 
us 3 TE - - 平均 相对 误差 (96) 
" 热流 密度 (W/m) 入 口 温度 (K) 热流 密度 (W/m) 入 口 温度 (K) 
1 3228.67 289.58 3300 290 
2 4266.56 294.54 4200 295 
3 5633.68 300.71 5700 300 1.069 
4 7367.75 305.67 7300 305 
5 9134.39 309.31 9200 310 


为 了 对 比 低 阶 反问 题 计算 模型 的 计算 效率 ， 本 ”实施 过 程 可 参见 文献 [15]。 表 4 给 出 了 共 斩 梯 度 法 的 
文 还 采用 了 经 典 的 共 轿 梯度 法 对 上 述 反 问题 进行 了 ”和 达 代 求解 结果 ， 其 中 收敛 精度 为 c=0.1%， 即 当前 过 
求解 。 由 于 本 文 讨论 的 重点 是 低 阶 反问 题 计算 模型 ， ” 代 层 次 的 结果 与 上 一 次 迭代 层次 的 结果 之 差 小 于 c 
内 此 ， 这 里 只 给 出 了 共 思 梯 度 法 的 计算 结果 ， 具 体 时， 结果 达到 收敛 。 


H 
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计算 所 采用 计算 机 的 CPU 频率 为 3.30 GHz， 内 
存 为 8.00 G. R 4 还 给 出 了 共 轿 梯度 法 与 低 阶 反问 
题 计 算 模型 计算 时 间 的 对 比 ， 由 表 4 可 明显 地 发 现 
低 阶 反问 题 模 型 的 优势 ， 尽 管 计算 结果 的 精度 有 较 
小 幅度 的 损失 ， 但 低 阶 反问 题 计 算 模型 相对 于 共 斩 
梯度 法 可 节省 1200 倍 以 上 的 计算 时 间 。 需 要 指出 的 


是 ， 这 里 的 计算 时 间 不 包括 样本 采集 及 POD Æ Kr 
求解 过 程 ， 尽 管 这 些 过 程 需 消耗 较 多 时 间 ， 但 对 于 
由 于 测量 到 的 数据 改变 而 需要 快速 计算 出 待 求 参 数 
时 ， 这 些 过 程 只 进行 一 次 ， 因 此 可 归结 为 前 期 准备 
工作 。 


表 4 RMR ARER 


Table 4 Computational results for conjugate gradient methodology 


温度 场 数 计算 结果 Jt pe HET 计算 时 间 (s) 
据 热流 密度 Wim) — Ar K) RAA — 低 阶 反问 题 计算 模型 共 斩 梯度 法 
1 3266.87 290.33 18 16.231 20160.263 
2 4229.66 294.16 17 15.489 19320.252 
3 5666.26 300.19 18 15.963 20230.446 
4 7272.98 305.21 19 16.660 21390.369 
5 9179.43 309.88 17 15.269 19080.298 
fluid flow with moving boundaries by proper orthogonal 
4 结 论 decomposition [J]. International Journal of Heat and Fluid 


本 文采 用 POD 技术 设计 了 一 般 工 程 反 问题 的 低 
阶 计 算 模型 ， 与 经 典 的 共 斩 梯 度 法 相 比 ， 低 阶 反 问 
题 计算 模 型 在 求解 过 程 中 不 涉及 原 物 理 问 题 的 控 表 
方程 ， 而 只 以 少量 的 POD 系数 作为 中 间 求 解 变量 ， 
在 大 幅度 简化 计算 过 程 的 同时 也 表明 了 该 方法 的 通 
用 性 。 以 圆 管内 强制 对 流 换 热 为 应 用 实例 ， 计 算 结 
果 表 明 ， 低 阶 反 问题 计算 模型 的 计算 结果 与 真实 结 
果 的 平均 相对 误差 约 为 1.0%， 计 算 速 度 比 经 典 的 共 
HERA 1200 倍 以 上 ， 证 明了 该 方法 的 正确 性 与 
可 行 性 ， 为 开发 工程 上 快速 预测 未 知 参数 的 技术 提 
供 了 理论 支撑 。 
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